Numerical implementation of a mixed finite element formulation for convection-diffusion problems

This document aims to the numerical solution of convection-diffusion problems in a fluid dynamics context by means of the Finite Element Method (FEM). It describes the classical finite element solution of convection-diffusion problems and presents the implementation and validation of a new formulation for improving the accuracy of the standard approach. On first place, the importance and need of numerical convection-diffusion models for Computational Fluid Dynamics (CFD) is emphasized, highlighting the similarities between the convection-diffusion equation and the governing equations of fluid dynamics for incompressible flow. The basic aspects of the finite element method needed for the standard solution of general convection-diffusion problems are then explained and applied to the steady state case. These include the weak formulation of the initial boundary value problem for the convection-diffusion equation and the posterior finite element spatial discretization of the weak form based on the Galerkin method. After their application to the steady transport equation a simple numerical test is performed to show that the standard Galerkin formulation is not stable in convection-dominated situations, and the need for stabilization is justified. Attention is then focused on the analysis of the truncation error provided by the Galerkin formulation, leading to the derivation of a classical stabilization technique based on the addition of artificial diffusion along the flow direction, the so-called streamline-upwind (SU) schemes. Next, a more general and modern stabilization approach known as the Sub-Grid-Scale (SGS) method is described, showing that SU schemes are a particular case of it. Taking into account all the concepts explained, a new mixed finite element formulation for convection-diffusion problems is presented. It has been proposed by Dr. Riccardo Rossi, a researcher from the International Center for Numerical Methods in Engineering (CIMNE), and consists on extending the original convection-diffusion equation to a system in mixed form in which both the unknown variable and its gradient are computed simultaneously, leading to an increase in the convergence rate of the solution. The formulation, which had not been tested before, is then implemented and validated by means of a multiphysics finite element software called \texttt{Kratos}. Eventually, the obtained results are analyzed, showing the improved performance of the mixed formulation in pure diffusion problems

Study of unstructured finite volume methods for the solution of the Euler equations

Development of an unstructured finite volume solver for the numerical solution of high-speed flows using the Euler equation set.This work deals with the numerical solution of inviscid compressible flows by means of the Euler equations. It focuses on the description of an unstructured finite volume method for these equations and its numerical application to solve external, two-dimensional steady problems. On first place, the standard formulation of the Euler equations is presented, reviewing the most important properties that characterize their mathematical behavior. The hyperbolic nature of the system is discussed, emphasizing the fundamental importance of taking into account the propagation of information in the flow field in order to obtain physically meaningful solutions, which also leads to a description of how the boundary conditions should be treated to avoid undesirable behaviors. To complete this presentation, a dimensionless form of the equations is derived, which provides substantial advantages to the numerical solution. The attention is then focused on the unstructured finite volume formulation, which is based on a central approximation of the fluxes at the volume interfaces. According to the need of properly accounting for the propagation of characteristic variables, the requirement to add artificial dissipation terms to the central discretization is justified. Then, two classical forms of artificial dissipation are defined, namely, the first-order upwind scheme and the Jameson-Schmidt-Turkel high-order model, detailing how to adapt the formulation of the dissipation terms to an unstructured mesh. Eventually, the time integration of the spatially discretized equations is assessed. With the objective of performing a practical implementation of the theoretical concepts studied, the development of a numerical solver is presented next, briefly describing the program structure and characteristics. After that, five different test cases are solved with the purpose of validating the code, consisting on two transonic flows around a NACA0012 airfoil and three supersonic examples, respectively around a NACA0012 airfoil, a double wedge airfoil and circular cylinder. The results obtained for each case are then analyzed and compared against reference solutions, showing an overall satisfactory performance of the solver developed

Anisotropic flow of charged hadrons, pions and (anti-)protons measured at high transverse momentum in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm NN}}=2.76 TeV

The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.Comment: 16 pages, 5 captioned figures, authors from page 11, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/186

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.Comment: 4 pages, 1 figur

Centrality dependence of charged particle production at large transverse momentum in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm{NN}}} = 2.76 TeV

The inclusive transverse momentum ($p_{\rm T}$) distributions of primary charged particles are measured in the pseudo-rapidity range $|\eta|<0.8$ as a function of event centrality in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV with ALICE at the LHC. The data are presented in the $p_{\rm T}$ range $0.15<p_{\rm T}<50$ GeV/$c$ for nine centrality intervals from 70-80% to 0-5%. The Pb-Pb spectra are presented in terms of the nuclear modification factor $R_{\rm{AA}}$ using a pp reference spectrum measured at the same collision energy. We observe that the suppression of high-$p_{\rm T}$ particles strongly depends on event centrality. In central collisions (0-5%) the yield is most suppressed with $R_{\rm{AA}}\approx0.13$ at $p_{\rm T}=6$-7 GeV/$c$. Above $p_{\rm T}=7$ GeV/$c$, there is a significant rise in the nuclear modification factor, which reaches $R_{\rm{AA}} \approx0.4$ for $p_{\rm T}>30$ GeV/$c$. In peripheral collisions (70-80%), the suppression is weaker with $R_{\rm{AA}} \approx 0.7$ almost independently of $p_{\rm T}$. The measured nuclear modification factors are compared to other measurements and model calculations.Comment: 17 pages, 4 captioned figures, 2 tables, authors from page 12, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/284

Charge separation relative to the reaction plane in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm NN}}= 2.76 TeV

Measurements of charge dependent azimuthal correlations with the ALICE detector at the LHC are reported for Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV. Two- and three-particle charge-dependent azimuthal correlations in the pseudo-rapidity range $|\eta| < 0.8$ are presented as a function of the collision centrality, particle separation in pseudo-rapidity, and transverse momentum. A clear signal compatible with a charge-dependent separation relative to the reaction plane is observed, which shows little or no collision energy dependence when compared to measurements at RHIC energies. This provides a new insight for understanding the nature of the charge dependent azimuthal correlations observed at RHIC and LHC energies.Comment: 12 pages, 3 captioned figures, authors from page 2 to 6, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/286

A note on comonotonicity and positivity of the control components of decoupled quadratic FBSDE

In this small note we are concerned with the solution of Forward-Backward Stochastic Differential Equations (FBSDE) with drivers that grow quadratically in the control component (quadratic growth FBSDE or qgFBSDE). The main theorem is a comparison result that allows comparing componentwise the signs of the control processes of two different qgFBSDE. As a byproduct one obtains conditions that allow establishing the positivity of the control process.Comment: accepted for publicatio

Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV

We report the first measurement of charged particle elliptic flow in Pb-Pb collisions at 2.76 TeV with the ALICE detector at the CERN Large Hadron Collider. The measurement is performed in the central pseudorapidity region (|$\eta$|<0.8) and transverse momentum range 0.2< $p_{\rm T}$< 5.0 GeV/$c$. The elliptic flow signal v$_2$, measured using the 4-particle correlation method, averaged over transverse momentum and pseudorapidity is 0.087 $\pm$ 0.002 (stat) $\pm$ 0.004 (syst) in the 40-50% centrality class. The differential elliptic flow v$_2(p_{\rm T})$ reaches a maximum of 0.2 near $p_{\rm T}$ = 3 GeV/$c$. Compared to RHIC Au-Au collisions at 200 GeV, the elliptic flow increases by about 30%. Some hydrodynamic model predictions which include viscous corrections are in agreement with the observed increase.Comment: 10 pages, 4 captioned figures, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/389

Measurement of charm production at central rapidity in proton-proton collisions at s=2.76\sqrt{s} = 2.76 TeV

The $p_{\rm T}$-differential production cross sections of the prompt (B feed-down subtracted) charmed mesons D$^0$, D$^+$, and D$^{*+}$ in the rapidity range $|y|<0.5$, and for transverse momentum $1< p_{\rm T} <12$ GeV/$c$, were measured in proton-proton collisions at $\sqrt{s} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. The analysis exploited the hadronic decays D$^0 \rightarrow$K$\pi$, D$^+ \rightarrow$K$\pi\pi$, D$^{*+} \rightarrow$D$^0\pi$, and their charge conjugates, and was performed on a $L_{\rm int} = 1.1$ nb$^{-1}$ event sample collected in 2011 with a minimum-bias trigger. The total charm production cross section at $\sqrt{s} = 2.76$ TeV and at 7 TeV was evaluated by extrapolating to the full phase space the $p_{\rm T}$-differential production cross sections at $\sqrt{s} = 2.76$ TeV and our previous measurements at $\sqrt{s} = 7$ TeV. The results were compared to existing measurements and to perturbative-QCD calculations. The fraction of cdbar D mesons produced in a vector state was also determined.Comment: 20 pages, 5 captioned figures, 4 tables, authors from page 15, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/307

Particle-yield modification in jet-like azimuthal di-hadron correlations in Pb-Pb collisions at sNN\sqrt{s_{\rm NN}} = 2.76 TeV

The yield of charged particles associated with high-$p_{\rm T}$ trigger particles ($8 < p_{\rm T} < 15$ GeV/$c$) is measured with the ALICE detector in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV relative to proton-proton collisions at the same energy. The conditional per-trigger yields are extracted from the narrow jet-like correlation peaks in azimuthal di-hadron correlations. In the 5% most central collisions, we observe that the yield of associated charged particles with transverse momenta $p_{\rm T}> 3$ GeV/$c$ on the away-side drops to about 60% of that observed in pp collisions, while on the near-side a moderate enhancement of 20-30% is found.Comment: 15 pages, 2 captioned figures, 1 table, authors from page 10, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/350